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Abstract 

The electron transfer in different solvents is investigated for systems consist- 
ing of donor, bridge and acceptor. It is assumed that vibrational relaxation 
is much faster than the electron transfer. Electron transfer rates and final 
populations of the acceptor state are calculated numerically and in an ap- 
proximate fashion analytically. In wide parameter regimes these solutions are 
in very good agreement. The theory is applied to the electron transfer in 
II2P — ZnP — Q with free-base porphyrin (H2P) being the donor, zinc por- 
phyrin (ZnP) the bridge, and quinone (Q) the acceptor. It is shown that the 
electron transfer rates can be controlled efficiently by changing the energy of 
the bridging level which can be done by changing the solvent. The effect of 

the solvent is determined for different models. 
PACS numbers: 31.70.Hq, 34.70.+e, 82.20.Rp 
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I. INTRODUCTION 



Electron transfer (ET) is a very important process in biology, chemistry and physics 
The most well known ET theory is the one of Marcus |0. Of special interest is the 
ET in configurations where a bridge (B) between donor (D) and acceptor (A) mediates the 
transfer. On this kind of ET we will focus in this paper. The primary step of ET in bacterial 
photosynthetic reaction centers is of this type |^ and a lot of work in this direction was done 
after the structure of the protein-pigment complex of the photosynthetic reaction centers 
of purple bacteria was clarified in 1984 . Many artificial systems especially self-organized 
porphyrin complexes have been developed to model this bacterial photosynthetic reaction 
center iijTl. 

Bridge-mediated ET reactions can occur via different mechanisms incoherent 



sequential transfer in which the bridge level is populated or coherent superexchange |14,15 



in which the mediating bridge level is not populated but nevertheless necessary for the 
transfer. Changing a building block of the complex |]9| JT6|JT7[| or changing the environment 
16| , p!8[| can modify which mechanism is mainly at work. Actually, there is still a discussion 



in literature whether sequential transfer and superexchange are limiting cases of one process 
p!9| or whether they are two processes which can coexist [0. To clarify which mechanism 
is present in an artificial system one can systematically vary the energetics of the complex. 
In experiments this is done by substituting parts of the complexes or by 



changing the polarity of the solvent [|1^]. Also the geometry and size of the bridge block can 
be varied, and in this way the length of the subsystem through which the electron has to be 



transfered P |ri| , |20| -p3| can be changed. 

Superexchange occurs due to coherent mixing of the three or more states of the system 
|l^,|T5|pl,|25| . The ET rate in this channel depends algebraically on the differences between 



the energy levels P,p!0| and decreases exponentially with increasing bridge length ||l4| , p3| , p5 



When incoherent effects such as dephasing dominate the transfer is mainly sequential |T2,23 



i. e., the levels are occupied mainly in sequential order [^|I2|,|T3|,|T6[ . The dependence on the 
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differences between the energy levels is exponential P,p!0|. An increase of the bridge length 
induces only a small reduction in the ET rate p,^^^^ . This is why sequential transfer 
is the desired process in molecular wires [p^ , p^ . 

In the superexchange case the dynamics is mainly Hamiltonian and can be described on 
the basis of the Schrodinger equation. The physically important results can be obtained by 
perturbation theory and, most successfully, by the semiclassical Marcus theory p. 

The complete system dynamics can be directly extracted by numerical diagonalization of 
the Hamiltonian p3| , |29|l . In case of sequential transfer the influence of an environment has to 
be taken into account. There are quite a few different ways how to include an environment 
modeled by a heat bath. The simplest phenomenological descriptions are based on the 
Einstein coefficients or on the imaginary terms in the Hamiltonian ||30|,|3T||, as well as on the 



Fokker-Planck or Langevin equations [0,^. The most accurate but also numerically most 
expensive way is the path integral method [^. This has been applied to bridge-mediated 
ET especially in the case of bacterial photosynthesis [^. Bridge- mediated ET has also 
been investigated using Redfield theory |]TB|,p3, by propagating a density matrix (DM) in 



Liouville space and other methods (e. g. P5| , PU| , |5^ -PB[]). 

The purpose of the present investigation is to present a simple, analytically solvable 
model based on the DM formalism [^,^ and apply it to a porphyrin-quinone complex 
which is taken as a model system for the bacterial photosynthetic reaction center. The 
master equation which governs the DM evolution as well as the appropriate relaxation 
coefficients can be derived from such basic informations as system-environment coupling 
strength and spectral density of the environment [p7| -|45[|. In the present model relaxation 
is introduced in a way similar to Redfield theory but in site representation not in eigenstate 
representation. A discussion of advantages and disadvantages these representations has been 
given elsewhere . The equations for the DM are the same as in the generalized stochastic 



Liouville equation (GSLE) model |^^^ for exciton transfer which is an extension of the 
Haken-Strobl-Reineker (HSR) model to a model with a quantum bath. Here we give 

an analytic solution to these equations. The present equations for the DM obtained are also 



similar to those of Ref. |TT| where relaxation is introduced in a phenomenological fashion but 
only a steady-state solution is found in contrast to the model introduced here. In addition 
the present model is applied to a concrete system. A comparison of the ET time with 
the bath correlation time allows us to regard three time intervals of system dynamics: the 
interval of memory effects, the dynamical interval, and the kinetic, long-time interval p9[| . 
In the framework of DM theory one can describe the ET dynamics in all three time intervals. 
However, often it is enough to find the solution in the kinetic interval for the explanation 
of experiments within the time resolution of most experimental setups, as has been done 



in Ref. |11,51|. The master equation is analytically solvable only for simple models, for 



example [^^. Most investigations are based on the numerical solution of this equation 
||T^ , |26| , ^ , |43[| . Here we perform numerical as well as approximate analytical calculations for 
a simple model. Since the solution can be easily obtained, the influence of all parameters 
on the ET can be examined. 

The paper is organized as follows. In the next section we introduce the model of a 
supermolecule which we use to describe ET processes. The properties of an isolated super- 
molecule are modeled in subsection 111 A|, as well as the static influence of the environment. 



The dynamical influence of bath fluctuations is discussed and modeled by a heat bath of 
harmonic oscillators in section |II B| . The reduced DM equation of motion (RDMEM) de- 
scribing the excited state dynamics is presented in subsection [II Q . In subsection flD\ the 



system parameter dependence on the solvent dielectric constant is discussed for different 
models of solute-solvent interaction. In subsection [II E| system parameters are determined. 
The methods and results of the numerical and analytical solutions of the RDMEM are pre- 
sented in section |T|. The dependencies of the ET rate and final acceptor population on the 
system parameters are given for the numerical and analytical solutions in subsection [IV A . 



The analysis of the physical processes in the system is also performed there. In subsection 
[IV B[ we discuss the dependence of the ET rate on the solvent dielectric constant for dif- 
ferent models of solute-solvent interaction and compare the calculated ET rates with the 
experimentally measured ones. The advantages and disadvantages of the presented method 



in comparison with the GSLE model ||47|j48[] and the method of Davis et ah are analyzed 
in subsection [IV C[ In the conclusions the achievements and possible extensions of this work 
are discussed. 



II. MODEL 

A. System Part of the Hamiltonian 

The photoinduced ET in supermolecules consisting of three sequentially connected molec- 
ular blocks, i. e., donor, bridge, and acceptor, (M = 1,2,3) is analyzed. The donor is not 
able to transfer its charge directly to acceptor because of their spatial separation. Donor 
and acceptor can exchange their charges only through B. In the present investigation the 
supermolecule consists of free-base tetraphenylporphyrin (112?) as donor, zinc substituted 
tetraphenylporphyrin (ZnP) as bridge, and p-benzoquinone as acceptor ||16|. In each of 
those molecular blocks we consider only two molecular orbitals (m = 0, 1), the highest oc- 
cupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) 
|53| . Each of these orbitals can be occupied by an electron or not, denoted by |1) or |0), 



respectively. This model allows us to describe four states of each molecular block, the neu- 
tral ground state |1)homo|0)lumo5 the neutral excited state |0)homo|1)lumo5 the positively 
charged ionic state |0)homo|0)lumo; and the negatively charged ionic state |1)homo|1)lumo- 
= |l)Mm(0|Mm, Cjum = |0)Mm(l|Afm, and = c\j^CMm dcscrlbc the creatlou, anni- 
hilation, and number of electrons in orbital Mm, respectively, while hm = J2m'^Mm gives 
the number of electrons in a molecular block. The number of particles in the whole super- 
molecule is conserved, Y^m^m = const. 

Each of the electronic states has its own vibrational substructure. As a rule for the 
porphyrin- containing systems the time of vibrational relaxation is two orders of magnitude 



faster than the characteristic ET time [16|. Because of this we assume that only the ground 



vibrational states play a role and we do not include the vibrational substructure. A com- 
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parison of the models with and without vibrational substructure has been given elsewhere 



Below we consider the evolution of single charge-transfer exciton states in the system. 
For the full description of the system one also should include photon modes to describe for 
example the fluorescence from the LUMO to the HOMO in each molecular block transferring 
an excitation to the electro-magnetic field. But the rates of fluorescence and recombination 
are small in comparison to other processes for porphyrin- type systems [^,^]. When fluores- 
cence does not have to be taken into account, all states except |D*BA) (M = 1), |D^B^A) 
(M = 2), and |D+BA^) (M = 3) remain essentially unoccupied, while those three take part 
in the intermolecular transport process. In this case the number of states coincides with the 
number of sites in the system and we label the states /i = 1,2,3 instead of {M, m}. 

For the description of the ET and other dynamical processes in the system placed in a 
dissipative environment we introduce the Hamiltonian 

H = Hs + Hb + Hsb, (1) 

where describes the supermolecule, Hb the dissipative bath, and Hsb their interaction. 
We are interested in the kinetic limit of the excited state dynamics here. For this limit 
we assume that the relaxation of the solvent takes only a very short time compared to the 
system times of interest. 

The effect of the solvent is twofold. On one hand the system dynamics is perturbed by 
the solvent state fluctuations, independent of the system states. Hsb shall only reflect the 
dynamical influence of the fluctuations leading to dissipative processes as discussed in the 



next subsection. On the other hand the system states are shifted in energy |56 



Hs = Ho + + V, (2) 

due to the static influence of the solvent which is determined by the relaxed value of the 
solvent polarization and in general also includes the non-electrostatic contributions such as 
van-der-Waals attraction, short-range repulsion, and hydrogen bonding |^^^. In Eq. 
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the energy of free and noninteracting blocks Hq = Y^Mm EMmnMm, is given by the energies 
Emtu of orbitals Mm in the independent electron approximation The Emtu are chosen 



to reproduce the ground-state-excited-state transitions e. g. D ^ D*, which change only a 
little for different solvents and are assumed to be constants here. To determine 



one starts from fully ionized double bonds in each molecular block [^], calculates the one- 
particle states and fills these orbitals with two electrons each starting from the lowest energy. 
By exciting, removing, adding the last electron to the model system one obtains the energy 
of the excited, oxidized, reduced molecular block in the independent particle approximation. 
The inter-block hopping term 

V = J:^,AV,1 + V,.) [in, - If + (n, - 1)' 

in Eq. (Q) includes the hopping operators V^i, = c^j^cmi, and the coherent couplings v^j,^. 
We assume f 13 = because there is no direct connection between donor and acceptor. The 



scahng of for different solvents is discussed in subsection [II D . 

The electrostatic interaction H^s scales like energies of a system of charges in a single 
or multiple cavity surrounded by a medium with static dielectric constant eg according to 
the classical reaction field theory [^]. Here we consider two models of scaling. In the first 
model each molecular block is in an individual cavity in the dielectric. For this case the 
electrostatic energy reads 

^cs = ^''(es)(^el + ^ion). (3) 

He\ = "^1^1,- l|e^ (47reor^)"^ 

takes the electron interaction into account while bringing an additional charge onto the 
block /i and thus describes the energy to create an isolated ion. This term depends on the 
characteristic radius of the molecular block. The interaction between the ions 
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depends on the distance between the molecular blocks r^^. Both distances and r^,^ are 
also used in Marcus theory H . The term H^i + //ion reflects the interaction of charges inside 



the supermolecule which is weakend by the reaction fleld according to the Born formula |61 



In the second model, considering the supermolecule as one object placed in a single cavity 
of constant radius one has to use the Onsager term |6l| . This term is state selective, it gives a 
contribution only for the states with nonzero dipole moment, i.e., charge separation. Deflning 
the static dipole moment operator as p = Y^{n^ — V){ny — V)ry,ye we obtain ifes = 5 p /ri3, 
with Onsager scaling 



B. Microscopic Motivation of the System-Bath Interaction and the Thermal Bath 

One can express the dynamic part of the system-bath interaction as 

HsB = -J d^rJ2 ■ A^(^)- (6) 



Here D^^{r) denotes the fleld of the electrostatic displacement at point r induced by the 
system transition dipole moment p^^ = p^iuiV^u + ^m*^) HII- The fleld of the environmental 
polarization is denoted as P{f) = J2n ^i'f^— ^n)dn, where dn is the nth dipole of the environ- 
ment and r„ its position. Only fluctuations of the environment polarization AP(r) influence 
the system dynamics. Averaged over the angular dependence the interaction reads |^ 



The dynamical influence of the solvent is described with a thermal bath model. The 



deviation A 



dr. 



of dn from its mean value is determined by temperature induced fluctuations. 



For unpolar solvents described by a set of harmonic oscillators the diagonalization of their 

8 



interaction yields a bath of harmonic oscillators with different frequencies uj\ and effective 
masses mx. In the case of a polar solvent the dipoles are interacting rotators as, e.g. used to 
describe magnetic phenomena |^2|,^. The elementary excitation of each frequency can again 



be characterized by an appropriate harmonic oscillator. So we use generalized coordinates of 



solvent harmonic oscillator modes Q\ = yh{2mxLUx) ^{ax + a^) for polar as well as unpolar 
solvents. The occupation of the ith state of the Ath oscillator is defined by the equilibrium 
DM px,ij = exp [~huJxi/{kBT)]6ij. 

All mutual orientations and distances of solvent molecules have equal probability. An 
average over all spatial configurations is performed. The interaction Hamiltonian (0) is 
written in a form which is bilinear in system and bath operators: 



SB 



J2Kxiat + ax) 



SsB (8) 



PfiuKx denotes the interaction intensity between the bath mode ax of frequency ux and 
the quantum transition between the LUMOs of molecules fi and u with frequency uj^^i, = 
{E^ — E^) /h. The scaling function Ssb refiects the properties of the solvent. Explicit 
expressions for the solvent infiuence are still under discussion in the literature p7| , |58| . 

C. Reduced Density Matrix Approach 

The interaction of the system with the bath of harmonic oscillators describes the ir- 
radiative energy transfer from the system to the solvent as modeled by Eq. (H). For the 
description of the dynamics we use the reduced DM which can be obtained from the full 



DM p by tracing over the environmental degrees of freedom a = TtbP [0 with the evolu- 



tion operator technique p!8| , |64| , restricting ourselves to the second order cumulant expansion 



Here we apply the Markov approximation, i.e., we restrict ourselves to the limit of long 
times. Furthermore, we replace the discrete set of bath modes with a continuous one. To 
do so one has to introduce the spectral density of bath modes J{uj) = nj^x K16{lj — ux)- 
Finally one obtains the following master equation 
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- {r^K [n{ujf,^,) + 1] + T^^n{uj^^) + T [n{uJ^,x) + 1] + rA;,n(u;A^)} (t«;a 

+ {Fa. [2n(cuA.) + 1] + T.a [2n{uj^x) + 1]} (Tak, (9) 

where n{ijj) = [exp {hu/kBT) — 1]^^ denotes the Bose-Einstein distribution. The damping 
constant 

r^. = sl^h~Miuj^,)pl^ (10) 

reflects the couphng of the transition — |z/) to a bath mode of the same frequency. It 
depends on the density of bath modes J at the transition frequency and on the tran- 
sition dipole moments Pf^^- A RDMEM of similar structure was used for the description 
of exciton transfer in the Haken, Strobl, and Reineker (HSR) model |]49|j50|| and the gen- 



eralized stochastic Liouville equation (GSLE) model |4^;^^. The HSR method originating 
from the stochastic bath model, is valid only in the high temperature limit The GSLE 



method P7| , ^ appeals to the quantum bath model with system-bath coupling of the form 
-^SB ~ V^V (jix + <^x) ) which modulates the system transition frequency. In Ref. p7|^ the 
equations for exciton motion are derived using the projection operator technique. Taking 
the different system-bath coupling we have derived the RDMEM which coincides with GSLE 



47| , |48[| . Both GSLE and our RDMEM are able to describe correctly finite temperatures. 
Below we neglect the last term of Eq. (|^) corresponding to the 7 term in the HSR and GSLE 
models because the rotating wave approximation (RWA) is applied. 

For the sake of convenience of analytical and numerical calculations we replace F^i, 
and the population of the corresponding bath mode n{u!^^) with the dissipative transitions 
dfiu = ^fj,u\n{u!fiu)\ and the corresponding dephasings 7^1, = J2k {d^^ + di^y) /2. With this one 
can express the RDMEM (Pf) in the form 



(11) 
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The parameters controlling the transitions between the selected states are discussed in sub- 
section p.1 E| . 



D. Scaling of the Damping Constants 

The relaxation coefficients Eq. (|1^) include the second power of the scaling function 
SsB because one constructs the relaxation term of Eq. (P) with the second power of the 
interaction Hamiltonian. The physical meaning of ifsB is similar to the interaction of the 
system dipole with a surrounding media. That is why it is reasonable to use the Onsager 
expression (|^) for S^b- In the work of Mataga, Kaifu, and Koizumi the interaction 
energy between the system dipole and the media scales in leading order as 



SB 



2(e, - 1) 2(eoo - 1) 



(13) 



2es + 1 2eoo + 1 

where eoo denotes the optical dielectric constant. From a recent paper of Georgievskii, Hsu, 
and Marcus [E7|l we extract F ~ — — for the multiple cavities model assuming e^j = eoo- 

L| |J eg Coo 



In terms of a scaling function it can be expressed as 

^SB = (l/es-l/eoo)^. (14) 



As we have already argued in 0] the coherent coupling v^^, between two electronic states 



scales with eg and e^o too, because a coherent transition in the system is accompanied by 
a transition of the environment state which is larger for solvents with larger polarity. As 
discussed above we neglect the vibrational substructure of each electronic state because the 
vibrational relaxation is about two orders of magnitude faster than the characteristic ET 
time. But in contrast to the model with vibrational substructure the present model does 
not involve any reaction coordinate. To reproduce the results of the more elaborate model 
with vibrational substructure one has to scale the electronic couplings Vf^u with the Franck- 
Condon overlap elements FFc(/i, 0, z/, 0) between the vibrational ground states of each pair 
of electronic surfaces 
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= vltyPFcifJ-, 0, 1/, 0) , (15) 

where f is the couphng of electronic states of the isolated molecule. For the calculation of 
the Franck-Condon factors one has to introduce the leading (mean) environment oscillator 
frequency ciJvib- Here c^vib = 1500 cm~^ is used which is similar to the frequency of the C-C 
stretching mode. With this scaling one implicitly introduces a reaction coordinate into the 
model. 



E. Model Parameters 

The dynamics of the system is controlled by the following parameters: energies of system 
states Ef^, coherent couplings f^j,, and damping constants T^^. 

On the basis of the spectral data and taking reference energy -Edba = we determine 
Eb*ba = 1-82 eV (in CH2CI2). We take the energy of the state with ET to Q from reference 
16| : -Ed+ba- = 1-42 eV [^]. Further Rempel et al. ||16| estimate the coupling of initially 



excited and charged bridge states (D*BA|iJ|D+B"A) = v^2 = 65 meV = 9.8 x 10^^ s~^ 
and the coupling of the two states with charge separation (D+B^A|iJ|D+BA^) = = 
2.2 meV = 3.3 x 10^^ s^^. The values of the couplings are essentially lower than the energy 
differences between the relevant system states 

huij > v^j. (16) 



This is the reason to remain in site representation instead of eigenstate representation p7 . 
The damping constants are found with help of the analytical solution derived at the end of 
the next section to be r2i = r23 = 2.25 x 10^^ s^^. The typical radius of the porphyrin ring 



is about r„ = 5 ± 1 A [|T^, while the distance r„,^ between the blocks of H2P — ZnP — Q 



reaches ri2 = 12.5 ± 1 A [^,0, r23 = 7 ± 1 A, ri3 = 14 ± 1 A. The main parameter which 
controls ET in a triad is the energy of the state -Ed+b-a- This state has a big dipole moment 
because of its charge separation and is therefore strongly influenced by the solvent. Because 
of the special importance of this value we calculate it for the different solvents as a matrix 
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element of the system Hamiltonian (2). The calculated values of the energies of the D"'"B A 
state for some solutions are shown in Table |. 

III. RESULTS 

The time evolution of the ET in the supermolecule is described by solving numerically 
and analytically Eqs. (|Tl])-(|T2|) with the initial condition of donor excitation with a tt pulse 
of appropriate frequency, i.e., the donor population is set to one. 

For the numerical simulation we express the system of Eqs. (|TT|)-(|T2|) in the form a = Aa, 
where ex is a vector of dimension 3^ for the model with 3 system states and the super-operator 
A is a matrix of dimension 3^ x 3^. We find an exponential growth of the acceptor population 

Psit) = P-sioo) [1 - exp i-kETt)] , (17) 

where for the solvent MTHF ^et ^ 3.59 x 10^ s"^ and P3(oo) ~ 0.9994. The population 
P2 which corresponds to charge localization on the bridge does not exceed 0.005. This 
means that in this case the superexchange mechanism dominates over the sequential transfer 
mechanism. Besides it ensures the validity of characterizing the system dynamics with ^3(00) 
and 

fcET = P3(00)^"[l-P3(t)]rft| \ (18) 

The alternative analytical approach is performed in the kinetic limit 

t > l/min(7^,). (19) 

In Laplace space the inequality ( [T9| ) reads s -C min(7^,^), where s denotes the Laplace 
variable. It is equivalent to replacing the factor l/{iuj^i, + 7^1, + s) in the Laplace transform 
of Eqs. (pn])-(p!2D with l/{iujfj,u + •y^u)- This trick allows to substitute the expressions (|12|) 
for non-diagonal elements of the DM into Eq. (|TT|). After this elimination we describe the 
coherent transitions to which the non-diagonal elements contribute by redefinition of the 
diagonal RDMEM (|TT|) 
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9^iy(^^iti + gutter (20) 



The transition coefficients g^i, contain dissipative and coherent contributions 

gi,v = d^^ + v^^v^^-i^u\h^ {ul^ + ^il^y^ \ (21) 

Now it is assumed that the bridge is not populated. This allows us to find the acceptor 
population in the form of Eq. (|T^) , where 



^ET = 5'23 + 92?,{9i2 - 932) (921 + fi'23) \ (22) 
P3(oo) = gi2923 [{921 + 923) kETV^ ■ (23) 

The value of Ti^,„ = Sl^h~'^J{uj^i^)p'j^^ can be found comparing the experimentally determined 
ET rate and Eq. (^21). To calculate J{ujp,v) would require a microscopic model. To avoid a 
microscopic consideration we simply take the same T for all transitions between excited 



states. The value of ET for H2P — ZnP — Q in MTHF is found by Rempel et al. |]T6| to 
be /cet = 3.6 ± .5 x 10^ s~^. If the bridge state has a rather high energy one can neglect 
thermally activated processes. V23 is negligibly small with respect to Vi2- In this case our 
result (I22I) reads 



k^i = ^^?2r2ir23 (^'^21 + r^i) ' (r2i + 123)"' . (24) 

With the relation r2i = r23 and the experimental /cet one obtains r2i = r23 — 2.25 x 
10^^ s^^. The fit of the numerical solution of Eqs. (0)-(12) to the experimental /cet in 



MTHF gives the same value. So the damping constants are fixed for a specific solvent and 
for other solvents they are calculated with the scaling functions. With this method the ET 
was found to occur with dominance of the superexchange mechanism with rates 4.6 x 10^ s~^ 
for CYCLO and 3.3 x 10^ s'^ for CH2CI2. 



IV. DISCUSSION 
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A. Sequential Versus Superexchange 



To discuss how the transfer mechanism depends on the change of parameters we calculate 
the system dynamics varying one parameter at a time. The dependencies of ^et and P3(oo) 
on Vi2, V23 and r2i, are shown in Figs. 2 and 3. The change of each parameter influences 
the transfer in a different way. 

In particular, k^T depends quadratically on f 12 from 10^^ s"^ to 10^^ s~^ in Fig. 2. Below 
it saturates at the lower bound /cet oc 3 x 10^ s~^. This corresponds to a crossover of the ET 
mechanism from superexchange to sequential transfer. But, due to the big energy difference 
between donor and bridge states the sequential transfer efficiency is extremely low. This 
is displayed by P3(oo) ~ 0. In the region vu ~ fis both mechanisms contribute to k^T- 
The decrease of P3(oo) in this region corresponds to coherent back transfer. The ET rate 
depends on ^23 in a similar way. At rather high values of f 12, V23 — 10^^ s~^ the relation (p!6D 
is no more valid. For this regime one has to use eigenstate instead of site representation 
because the wavefunctions are no more localized . 



The variation r2i, r23 near the experimental values shows similar behavior of A;ET(r2i) 
and /cET(r23) (see Fig. 3). Here we independently vary r2i and r23. Both, A;ET(r2i) and 
^ET(r23) increase linear until the saturation value 7 x 10^ s'^ at F > 10^^ s~^ is reached. 
There is qualitative agreement between the numerical and analytical values. In Eq. (|TBp 
infinite time is approximated by 10^^ s and so one cannot obtain ET rates lower than this 
limit. 

The physical meaning of the ET rate dependence on F seems to be transparent. At 
small values of F a part of the population coherently oscillates back and forth between the 
states. The increase of the dephasing 7^1, quenches the coherence and makes the transfer 
irreversible. So transfer becomes faster up to a maximal value. For the whole range of F, 
depopulations ^21, ^23 and thermally activated transitions du, ^32 always remain smaller 
than the coherent couplings, therefore they do not play an essential role. 

Next, the similarity of the dependencies on F21 and F23 will be discussed on the basis 
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of Eq. (^). In the limit kBT/hu^u — > thermally activated processes with u^u < vanish 
and so |n(ti;^,^)| = 0, while depopulations with u;^,^ > remain constant |n(ti;^,^)| = 1. 
The condition cj^j^ ^ 7^j, allows us to neglect 7^^ in comparison with cj^^,. With these 
simplifications Eq. (^) becomes 

^ET — r2ir23 (r21 + r23)"' {v!,/ul + Vjjuji,) , (25) 

i.e. symmetric with respect to and r23. 

To the largest extent the mechanism of transfer depends on the bridge energy Ed+b-a 
as presented in Fig. 4. In different regions one observes different types of dynamics. For 
large bridge energy £'21 = -^d+b-a ~ -E'd*ba ^ the numerical and analytical solutions do 
not differ from each other. The transfer occurs with the superexchange mechanism. The 
ET rate reaches a maximal value of 10^^ s^^ for low bridge energies. 

While the bridge energy approaches the donor energy the sequential transfer starts to 
contribute to the ET process. The traditional scheme of sequential transfer is obtained when 
donor, bridge, and acceptor levels are arranged in a cascade. In this region the analytical 
solution need not coincide with the numerical solution because the used approximations are 
no more valid. For equal bridge and acceptor energies /cet displays a small resonance peak 
in Fig. 4(a). When the bridge energy is lower than the acceptor energy the population gets 
trapped at the bridge. The finite /cet for £'21 < does not mean ET because ^3(00) 0. 
For the dynamic time interval t < 7~J a part of the population tunnels force and back to 
the acceptor with k-^T- The analytical solution ( p2D gives a constant rate for the regime 
E21 < E31, while the numerical solution of Eqs. (|lT])-(|12D is instable, because such coherent 
oscillations of population cannot be described by Eq. ([I7|) and /cet cannot be fitted with 
Eq. ([T8|). In Fig. 4 the regime E21 < -E31 occurs for small E21 while E^i is kept constant and 
for large E^i while E21 remains constant. 

The energy dependence of the final population has a transparent physical meaning for the 
whole range of energy. A large bridge energy ensures the transition of the whole population 
to the acceptor. In the intermediate case, when the bridge has the same energy as the 
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acceptor, final population spreads itself over these two states ^3(00) = .5. Lowering the 
bridge even more the whole population remains on the bridge as the lowest state of the 
system. The dependence of the ET rate on the acceptor energy E^i = -Ed+ba- ~ -£^d*ba 
in Fig. 4 remains constant while the acceptor energy lies below the bridge energy. Increase 
of E31 up to E21 = 1.36 eV gives the maximal k-^T oc r2i. When E^i increases further 
the acceptor becomes the highest level in the system and therefore the population cannot 
remain on it. 



B. Different Solvents 

For the application of the results to various solvents and comparison with experiment one 
should use the scaling for energy, coherent couplings, and damping constants as discussed 
above. The combinations of the energy scaling in subsection [11 A| and damping constants 
scalings in subsection |11D| are represented in Fig. 5. An increase in eg from 2 to 4 leads 



to an increase of the ET rate, no matter which scaling is used. Further increase of eg 
induces saturation for the Onsager-Mataga scaling and even a small decrease Within the 
applied approximations an increase in the solvent polarizability and, consequently, of its 
dielectric constant lowers the bridge and acceptor energies and increases the system-bath 
interaction and, consequently, the relaxation coefficients. It induces a smooth rise of the ET 
rate for the Onsager-Mataga scaling. On the other hand large eg leads to essentially different 
polarisational states of the environment for the supermolecule states with different dipole 
moment. This reduces the coherent couplings, see Eq. (pIS]), leading for the Born-Marcus 
scaling to a small decrease of /cet for large eg. The ET rate with this scaling comes closer to 
the experimental value fcET(e^^^*"^^)- This gives a hint that the model of individual cavities 
for each molecular block is closer to reality than the model with a single cavity for the whole 
supermolecule. 

Below we consider Born scaling Eq. (^) for the system energies and Marcus scaling 
Eq. (JT^) for the damping constant to compare the calculated ET rates with the measured 
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ones. For the solvents CYCLO, MTHF, and CH2CI2 one obtains the relative bridge energies 
E21 = 1.77 eV, 1.36 eV, and 1.30 eV, respectively. 

The calculated ET rate coincides with the experimental value ||T6[ for H2P — ZnP — Q in 
CYCLO, see table |I[ For CH2CI2 the numerical ET rate is approximately 30% faster than 
the experimental value. It has to be noted that a value for the damping rates can be chosen 
such that the calculated curve almost passes through all three experimental error bars. On 
the other hand an error in the present calculation could be due to (i) absence of vibrational 
substructure of the electronic states in the present model; (ii) incorrect dependence of system 
states energies on the solvent properties; (iii) opening of additional transfer channels not 
mentioned in the scheme shown in Fig. 1. Each of these possibilities needs some comments. 

ad (i): The incorporation of the vibrational substructure will result in a complication 
of the model [|13|,^ giving a more complicated ET rate dependence on the energy of the 
electronic states and dielectric constant. It should yield the maximal ET rate for nonequal 
energies of electronic states, namely for the activationless case when the energy difference 
equals the reorganization energy. For a comparison of the models with and without vibra- 
tional substructure see Ref. 

ad (ii): Effects as the solvation shell [^] do need a molecular dynamics simulation. 
The total influence of the solvent is, probably, reflected in an energy shift between the 
spectroscopically observable states -Ed*ba and -Edb*a [0- 

ad (iii): A solvent with large eg can bring high- lying system states closer to the ones 
included in Fig. 1. E.g., because of its larger dipole moment, |D^B+A) is strongly influenced 
by the solvent. 



C. Comparison with similar theories 

As discussed above the RDMEM are very similar to those of the GSLE model. This is 
an extension of the HSR theory in which a classical bath is used and for which analytical 
solutions are available P5| , ^ . We are not aware of any analytical solution of the GSLE 
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model as presented here. Also this model has not been applied to similar ET processes. 

The numerical steady-state method used by Davis et al. ||11|] is an attractive one due 
to its simplicity, but unlike our method it is not able to give information about the time 
evolution of the system. We use a similar approach derived within a Redfield-like theory. 
But we consider dephasing and depopulation between each pair of levels. In contrast Davis 
et al. incorporate relaxation phenomenologically only to selected levels, dephasing 7 occurs 
between excited levels, while depopulation k takes place only for the sink from acceptor to the 
ground state. The advantage of the approach of Davis et al. is the possibility to investigate 
the ET rate dependence for the bridge consisting of more than one molecular block. This 
was not the goal of the present work but it can be extended into this direction. We are 
interested in the ET in a concrete molecular complex with realistic parameter values and 
realistic possibilities to modify those parameter. Our results as well as the results of Davis 
show that ET can occur as coherent (with the superexchange mechanism) or dissipative 
process (with the sequential transfer mechanism). 

V. CONCLUSIONS 

We have performed a study of the ET in the supermolecular complex H2P — ZnP — 
Q within the DM formalism. The determined analytical and numerical ET rates are in 
reasonable correspondence with the experimental data. The superexchange mechanism of 
ET dominates over the sequential transfer. We have investigated the stability of the model 
varying one parameter at a time. The qualitative character of the transfer is stable with 
respect to a local change of system parameter. The crossover between the two transfer 
mechanisms can be induced by lowering the bridge energy. The relation of the theory 
presented here to other theoretical approaches to ET has been discussed. 

The calculations performed in the framework of the present formalism can be extended 
in the following directions: (i) Considerations beyond the kinetic limit. The vibrational 
substructure has to be included into the model as well as solvent dynamics and, probably, 
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non-Markovian RDMEM. (ii) Enlargement of the number of molecular blocks in the complex, 
(iii) Initial excitation of states with rather high energy should open additional transfer 
channels. 
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FIGURES 

FIG. 1. Schematic presentation of the energy levels in the H2P — ZnP — Q complex. The 
three states in the boxes play the main role in FT which can happen either sequentially or by a 
superexchange mechanism. Dashed lines refer to sequential transfer, curved solid line to superex- 
change, dot-dashed to energy transfer followed by FT, dotted line optical excitation, and straight 
solid lines either fluorescence or irradiative recombinations. 

FIG. 2. Dependence of the FT rate (a) and the final acceptor population (b) on the co- 
herent couplings V12 (triangles and dashed line, V23 = = 2.2 meV), V23 (dots and solid line 
^12 = "^12 = 65 meV). Symbols represent numerical solution of Fqs. (p!l|)-(|T^), lines analytical 



solution (22)- (PI). 



FIG. 3. Dependence of the FT rate on the damping constants (triangles and dashed line), 
r23 (dots and solid line). The other parameters correspond to H2P — ZnP — Q in MTHF. Symbols 
represent numerical and lines analytical solution. 

FIG. 4. Dependence of the FT rate (a) and the final acceptor population (b) on the energy of B 
E21 = -Ed+b-a — -E'd'BA (triangles and dashed line, E31 = —0.4 eV) and A E31 = E'^+g^- — £'d*ba 
(dots and solid line, £'21 = 1-36 eV). Symbols represent numerical and lines analytical solution. 
vi2 = 65 meV, ^23 = 2.2 meV, T21 = T23 = 2.25 x lO^^ 3"^ 

FIG. 5. The FT rate /cet versus dielectric constant. The energies of the bridge and acceptor 
scale in accordance with Born expression (^ (triangles and solid line), Onsager expression 
(diamonds and dashed line). Coherent couplings and damping constants scale in accordance with 
Mataga's expression (^) (diamonds and dashed line), Georgievskii-Marcus expression (14) (tri- 



angles and solid line). Symbols represent numerical and lines analytical solution. Solid crosses 
with error bars give experimental values Note that by using a different value of the damping 
parameter F curves can be calculated which almost pass through all three experimental error bars. 
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TABLES 

TABLE L Energy of the charged bridge state in different solvents and corresponding ET rates 
(for calculations Born scaling @ of energy and Marcus scaling (|T^ ) of dissipation are used) . MTHF 
denotes 2-methyl-tetrahydrofuran, and CYCLO denotes cyclohexane. 

Solution CH2CI2 MTHF CYCLO 



es, |63] 9.08 6.24 2.02 

^D+B-A, eV 3.12 3.18 3.59 

/CET, lO'^ s~\ num. 33 36 0.46 

kET, 10^ s-\ an. 33 36 0.46 

/cET, 10^ s-\ exp. m 23 ±5 36 ±5 + 3 



28 



— Born-Marcus, an. 

— Onsager-Mataga, an. 
X Experiment 

A Born-Marcus, num. 
♦ Onsager-Mataga, num. 



